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Abstract. A massive long-lived r neutrino in the MeV regime modifies the primordial light-element 
abundances predicted by big-bang nucleosynthesis (BBN) calculations. This effect has been used to 
derive limits on m„ r . Because recently the observational situation has become somewhat confusing 
and, in part, intrinsically inconsistent, we reconsider the BBN limits on m VT . To this end we use our 
newly developed BBN code to calculate the primordial abundances as a function of m VT and of the 
cosmic baryon density r\. We derive concordance regions in the r\-m VT -plane for several sets of alleged 
primordial abundances. In some cases a concordance region exists only for a nonvanishing to„ t . At the 
present time BBN does not provide clear evidence for or against a r neutrino mass. 
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1 Introduction 

The theory of big-bang nucleosynthesis (BBN) describes 
the creation of the light elements hydrogen, helium, 
lithium, and beryllium in the early universe. The stan- 
dard version of this theory is based on three main as- 
sumptions. 1. The evolution of the universe is described 
by the homogeneous and isotropic Friedmann- Robertson- 
Walker model. 2. The standard model of particle physics 
is used. 3. The lepton asymmetry of the universe is of the 
same order as the baryon asymmetry so that the chemical 
potential of the neutrinos in the early universe is small 
relative to the temperature. 
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Fig. 1. Predicted values for the primordial abundances of the 
light elements assuming 2 (dashed), 3 (solid), and 4 (dotted) 
light neutrino families. The boxes indicate observationally in- 
ferred limits (Sect. 2). 



Because the nuclear reaction rates and the neutron 
half-life are well known, there remains only one free pa- 
rameter, namely the present-day cosmic number density 
r\ of baryons relative to photons. For a given r\ it is then 
possible to solve numerically a nuclear network in the ex- 
panding universe which gives the primordial light element 
abundances. In Fig. 1 we show the results of such a cal- 
culation for 10~ 10 < rj < 10~ 9 , based on our newly de- 
veloped BBN-code which is documented in Rehm (1996). 



The results agree very well with those of other groups 
(e.g. Yang et al. 1984; Smith et al. 1993). 

In the past, BBN has become a cornerstone of the 
hot big bang scenario in that one could find a concor- 
dance interval for r/ where all observed light-element abun- 
dances agreed with the calculated values. Further, BBN 
has been used to constrain new particle-physics mod- 
els, notably the properties of neutrinos, and thus has 
become a powerful tool in the field of particle astro- 
physics (e.g. Malaney & Mathews 1993). In particular, if 
m„ T , the mass of the r neutrino, is larger than a few 
100 kcV the cosmic expansion rate during the BBN epoch 
is modified enough to cause significant changes in the pre- 
dicted light-element abundances. Therefore, one can de- 
rive constraints on a putative mass for the r neutrino 
or other neutral leptons that might be present in the 
early universe (Dicus et al. 1978, Miyama & Sato 1978, 
Kolb & Scherrer 1982, Dolgov & Kirilova 1988). 

Since this earlier works, the m VT question has been 
revisited by many groups 

of authors (Kolb et al. 1991, Dolgov & Rothstein 1993, 
Kawasaki et al. 1994, Dolgov et al. 1996, Hannestad & 
Madsen 1996a,b, Fields et al. 1997). The measured width 
of the Z° boson shows that there are exactly three neu- 
trinos which are lighter than half the Z° mass so that the 
post-1990 works indeed focus on the question of a r neu- 
trino mass rather than on the number of light neutrino 
flavors or the mass of some arbitrary new lepton. One 
motivation for the post-1990 works was to understand the 
impact of certain approximations for the kinetic treatment 
of massive neutrino freeze-out. Another motivation was to 
include the latest observationally inferred primordial light- 
element abundances in the derivation of the constraints. 
While the treatments by the various groups differ in de- 
tail, they all agree that m Vr above a few 100 keV is for- 
bidden. A mass beyond a few 10 MeV would again be al- 
lowed by BBN, but is experimentally ruled out (ALEPH- 
collaboration Buskulic et al. 1995, Passalacqua 1996). 

Meanwhile, the observational situation has changed in 
that the first direct measurements of deuterium in quasar 
absorption systems have become available. Unfortunately, 
the results by different authors and from different sys- 
tems do not agree with each other, and do not necessarily 
agree with the classic interstellar-medium determination 
of the primordial abundance. In addition, a new determi- 
nation of the primordial 4 He abundance yields a signifi- 
cantly higher value than the previous "standard" result. 
Therefore, at the present time it is not clear which set of 
observational data (if any of the currently available ones) 
should be used to compare with the BBN calculations. It 
is not even clear whether there indeed exists any concor- 
dance interval for 77, and if so, where exactly it lies. This 
constitutes the current debate about a possible "BBN cri- 
sis" (see Hata et al. 1995, Copi et al. 1995). 

Therefore, while the most recent papers on the m Vr 
question have focused on various aspects of the kinetic 
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Table 1. Observationally inferred primordial light-element abundances. "Measured" values and the errors are taken from the 
given papers in which references to the original observations can be found. The errors are systematic except for those with the 
superscript "stat" (for statistical). For 7 Li two sets of systematic errors are given as explained in the text. 



Abundance Measurement Adopted Range Label Reference 

% 0.234 ± 0.003 stat ± 0.005 0.223-0.245 FHe Olive & Scully (1996) 

0.243 ± 0.003 stat 0.237-0.249 h 4 He Izotov et al. (1997) 

D/H x 10 5 1.6 ± 0.2 1.8-4.4 ISM Linsky et al. (1993), Dearborn et al. (1996) 

19 ±4 15-23 hQAS Rugers & Hogan (1996) 

2.3 ± 0.3 stat ± 0.3 1.4-3.2 1QAS Tytler et al. (1996) 

7 Li/HxlQ 10 1.6±0.07 stat l°, 4 l^ 0.9-3.7 7 Li Fields et al. (1996) 



treatment for these particles during and after freeze-out, 
the main current source of uncertainty for m„ T limits is 
actually the observational situation. The main purpose of 
the present note is to reconsider BBN limits on m„ T in 
the light of the current debate about the correct primor- 
dial light-clement abundances. In this regard our study is 
motivated by a similar discussion of the allowed number 
of massless neutrino species (Cardall & Fuller 1996). 

These authors have calculated the primordial light- 
element abundances as a function of the cosmic baryon-to- 
photon ratio 77, and of the equivalent number of massless 
neutrino species N e g which characterizes the energy den- 
sity and thus the expansion rate of the universe at the 
time of BBN. They have then derived detailed concor- 
dance regions in the 77-./V e fT-plane where different sets of 
observational data agree with the calculated abundances. 
However, these results cannot be trivially translated into 
concordance regions in the r\-m Vr -plane because the con- 
tribution of a massive tau neutrino to the expansion rate 
cannot be expressed in terms of an equivalent N g, as we 
will discuss below. Therefore, it is best to construct such 
concordance regions directly in the r\-m VT -plane. Before 
the onset of the current observational debate this sort of 
analysis was performed by Kawasaki et al. (1994). 

A rather complete numerical treatment of the v T kinet- 
ics was performed in the works of Kawasaki et al. (1994) 
and Hannestad & Madsen (1996a, b). Together with the 
recent studies by 

Dolgov et al. (1996) and Fields et al. (1997) they indicate 
that the final concordance regions in the r\-m Vr -plane arc 
rather insensitive to fine points of the v T kinetics. There- 
fore, we limit ourselves to the simplest approximation in 
which kinetic equilibrium is assumed for v T throughout so 
that an analytically integrated Boltzmann collision equa- 
tion can be used. With this approach we can reproduce 
previous results with sufficient accuracy. 

Of course, on a timcscale much longer than the BBN 
epoch (approximately 200 seconds), the r neutrinos have 
to decay into relativistic particles to prevent the uni- 
verse from becoming overdominatcd by r neutrinos to- 
day (e.g. Kolb & Turner 1990). Wc will not discuss the 
impact of neutrino decays on BBN. Extensive previous 



discussions include those by Kawasaki et al. (1994) and 
Dodelson et al. (1994). 

We begin our study in Sect. 2 with a short review of 
the current status of the observations. In Sect. 3 we de- 
scribe the influence of a massive r neutrino on the cosmic 
expansion rate, which is calculated by solving a simplified 
Boltzmann equation. In Sect. 4 the influence of a massive 
t neutrino on the primordial abundances is studied. We 
derive concordance regions in the r\-m VT -plane for different 
sets of observational data. Finally, we discuss and summa- 
rize our results in Sect. 5. Some details about the neutrino 
reaction rates are given in Appendix A. 

2 Observations 

In order to derive concordance intervals for the cosmic 
baryon density, or in order to derive BBN limits on novel 
particle-physics models, one must compare the calculated 
light-element abundances with observations. For most el- 
ements the difficulty is to derive the primordial abun- 
dances from values which are measured today in our vicin- 
ity. Even for He, where this appears to be relatively 
straightforward and model- independent, some controversy 
has emerged. Therefore, we begin with a discussion of the 
current observational situation. 

Hclium-4 

A large number of measurements exist for the 4 He as 
well as for the O and N abundances in low-metallicity 
HII regions. Using O or N as tracers, the 4 He abun- 
dance is extrapolated to zero metallicity to obtain 
the primordial value. We quote the usual "standard" 
value ("low 4 He" labeled l 4 He in Table 1) according to 
Olive & Scully (1996). Their results are very similar to 
those of Olive & Steigman (1995); the values and errors 
differ at most by 2 x 10~ 3 . In contrast, Izotov et al. (1997) 
have found a significantly higher value (h 4 He in Table 1) 
using a new set of emissivities and collisional enhancement 
factors for the analysis of the observed spectral lines and 
by using other selection criteria for the HII regions. Their 
method and interpretation of the differences, in turn, has 
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been critiqued by Olive et al. (1997), who reconfirm the 
above "low" value. 

We emphasize that we do not give more weight to any 
of the two values; we merely quote what is currently dis- 
cussed in the literature. In view of the disagreement be- 
tween the lower deuterium value inferred from quasar ob- 
servations (see below) it has been proposed that the sys- 
tematic uncertainties might be substantially larger than 
commonly assumed (Sasselov & Goldwirth 1995). They 
argue that the primordial mass fraction of 4 He might be 
as high as 0.26. 

Deuterium 

For deuterium, there is the canonical value from observa- 
tions of deuterium in the solar system (Geiss 1993) and 
the interstellar medium (Linsky et al. 1993, 
Linsky & Wood 1996, Piskunov et al. 1997) which pro- 
vide a lower limit on the primordial deuterium abun- 
dance. Based on the theory of galacto-chemical evolu- 
tion for deuterium, an upper limit — usually about a fac- 
tor of 2-3 higher (cf. Edmunds 1994) — can be inferred. 
The most recent analysis is that of Dearborn et al. (1996) 
which is labeled "ISM" in Table 1. In all canonical mod- 
els of galacto-chemical evolution a strong overproduction 
of 3 He with respect to observations in the solar sys- 
tem and the ISM is found. As long as this problem re- 
mains unsolved, limits based on an analysis of the sum 
of deuterium and 3 He should not be used to constrain 
the primordial abundances. The value quoted here re- 
lies on deuterium depletion only. Non-canonical models 
for stellar evolution, including additional mixing on the 
Red Giant Branch avoid a net production of 3 He in low 
mass stars and might provide a solution to this problem 
(Charbonnel 1995, Weiss et al. 1996). 

In addition, deuterium is the only case where one be- 
lieves to observe a primordial abundance directly. Deu- 
terium lines were found in extragalactic HII clouds which 
lie on the line of sight to high rcdshift quasars (Quasar 
Absorption Systems or QAS). Unfortunately, the val- 
ues found in different systems by different authors do 
not agree. Tytler et al. (1996) argue that the correct 
value is one order of magnitude smaller ("low QAS 
value" labeled 1QAS in Table 1) than the one originally 
found by Songaila et al. (1994) and Carswell et al. (1994). 
Their high value (hQAS in Table 1) was recently con- 
firmed in an analysis by Rugcrs & Hogan (1996). The high 
QAS deuterium abundance would allow for an r\ concor- 
dance interval at relatively low values when combined 
with the lower 4 He abundance, but poses severe prob- 
lems for the theory of galacto-chemical evolution of deu- 
terium and 3 He (e.g. Steigman 1994, Galli et al. 1995). 
Even the non-canonical models of Charbonnel (1995) and 
Weiss et al. (1996) cannot make the high deuterium value 
consistent with the present-day 3 He observations. It is ar- 
gued by Songaila et al. (1994), Carswell et al. (1994), and 



Tytler et al. (1996) that the high QAS deuterium value 
is in fact not deuterium, but another hydrogen cloud 
at a different position which mimics the deuterium line, 
and that the systems used are not sensitive for abun- 
dances as low as 1CP 5 . This interpretation is rejected 
by Rugers & Hogan (1996). Clearly, more high rcdshift 
quasar spectra are needed to reveal the true primordial 
deuterium abundance. 

Lithium-7 

The chemical evolution of lithium in stars is not well un- 
derstood and therefore the derivation of the primordial 
value is much more uncertain. We quote here a recent 
value from Fields et al. (1996), based on the observations 
by Molaro et al. (1995), for the so-called "Spite plateau", 
which is the minimal 7 Li abundance found in old Pop II 
halo stars of our galaxy. Since it is nearly constant in 
stars with a surface temperature of more than 5500 K 
and metallicities lower than about 5 % solar, it is usu- 
ally argued that this is the primordial value. Because the 
abundance depends on the modeling of the stellar atmo- 
sphere, the error is dominated by systematic effects (first 
set of systematic errors in Table 1). Furthermore, it has 
been suggested that 7 Li might have been depleted signif- 
icantly (Deliyannis et al. 1990, Pinsonneault et al. 1992), 
but the detection of the more fragile isotope 6 Li in two of 
these stars may be taken as an indication against strong 
depletion (Steigman et al. 1993). The possible depletion 
effect is represented by the second set of systematic errors 
in the 7 Li abundance. 

The current observational evidence is summarized in Ta- 
ble 1. For each observation we also give a range for the 
relevant abundance which we have adopted to compare 
with the calculations below. The adopted range has been 
derived by adding the stated systematic error linearly with 
the 2 a statistical one. For lithium, both systematic errors 
have been linearly added. These procedures are somewhat 
arbitrary, but represent the common practice in this field. 

3 Massive long-lived r neutrinos in the early uni- 
verse 

As a further ingredient for our study of possible BBN lim- 
its on m„ T we need to specify the contribution of a massive 
neutrino to the expansion rate of the early universe. In the 
radiation dominated phase the expansion rate is governed 
by the Friedmann equation without a curvature term 1 



1 We always use the system of natural units where h — c = 
k B = 1. 
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Fig. 2. Evolution of the parameter g*(t) for a Majorana r 
neutrino with m Vr = 10, 20, and 100 MeV (solid lines). For 
comparison, the evolution of g t for 2 and 3 massless neutrino 
families is also shown (dotted lines). 



Here, H = R/R is the expansion rate, mpi the Planck 
mass, and g the energy density in radiation. It is com- 
monly expressed in the form 



two neutrinos are taken to be massless. Thermodynamic 
equilibrium is maintained as long as reactions such as 

v T + v T < — ► a + a (4) 

are fast compared to the expansion of the universe. Here, 
a stands for any fermion which is kincmatically allowed. 
As the universe expands, this reaction slows down and 
eventually the number density of r neutrinos freezes out 
and stays constant thereafter (Fig. 3). 
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where T is the cosmic temperature and g* the effective 
number of thermally excited relativistic degrees of free- 
dom. This parameter is a function of temperature and 
thus of cosmic time. Just before the BBN epoch one finds 
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I {g e + 2N V ) = 10.75 



(3) 



with g-y = 2 for the photons, g e = 4 for the electrons and 
positrons, and N u = 3 for the number of neutrino families 
which are assumed to be massless. 

A particle species no longer contributes when T falls 
below its mass because then its number density is sup- 
pressed and therefore the contribution to the cosmic en- 
ergy density in the early radiation dominated phase of the 
universe is negligible. The evolution of in the temper- 
ature range 10 2 to 10~ 2 MeV is shown in Fig. 2 (dashed 
lines) for a scenario with 2 and 3 light neutrino families. 
At a temperature of about 1 MeV the electrons get nonrel- 
ativistic and annihilate to photons. Therefore they drop 
out of Eq. (3), g* decreases, and their entropy is trans- 
ferred to the photons, which dominate the energy density 
from now on. 

Massive neutrinos contribute differently to the cosmic 
energy density than massless ones. Therefore, the expan- 
sion rate as a function of temperature and thus the time- 
temperature relation is modified, leading to changed light- 
element abundances in a BBN calculation. We examine a 
scenario where the r neutrino is massive whereas the other 
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Fig. 3. Number density (upper panel) and energy density 
(lower panel) as a function of the cosmic temperature for a 
10 MeV Majorana r neutrino. The contribution of a massless 
neutrino is shown for comparison. The densities are relative to 
those of photons and normalized to 1 at high temperatures. At 
e + e~ annihilation the photon density increases and the rela- 
tive neutrino densities take a downward step as indicated. The 
curves marked "Equilibrium" represent the v T densities if the 
v T v T annihilations would not freeze out. 

The process of freeze-out is described by the Boltz- 
mann collision equation. This intcgro-differential equation 
is very difficult to solve in general, but it can be simplified 
by making the following approximations: 1. The annihila- 
tion products arc in thermodynamic equilibrium due to 
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effective interactions with the cosmic plasma. 2. Classical 
instead of quantum statistics is used, i.e. a Maxwcll- 
Boltzmann distribution instead of a Fcrmi-Dirac or a 
Bosc-Einstein one. 3. Even after annihilations the t neu- 
trinos are kept in kinetic equilibrium by effective elastic 
scatterings with electrons and massless neutrinos which 
are far more abundant. Other authors have dropped some 
or all of these assumptions, but the results do not change 
drastically, as we will discuss in Sect. 4. Because of the 
assumed kinetic equilibrium, the distribution function of 
the massive t neutrinos can be expressed in terms of the 
so-called pseudo chemical potential z%, 

Ur = CX P - Y) ■ ^ 

Furthermore, the Boltzmann equation is integrated over 
the phase space of one of the incoming r neutrinos. Be- 
cause we assumed classical statistics, this can be done in 
a closed form and one arrives at one single ordinary dif- 
ferential equation for the v T number density, 

-g. + Wn VT = (av) [n„ r (t) 2 - n£ (t) 2 } . (6) 

Here, H is the expansion rate, n„ T the actual and n e °^ the 
equilibrium number density of r neutrinos. The thermally 
averaged reaction rate (av) is given in Appendix A. 

For a Dirac neutrino we are confronted with another 
difficulty because it has four degrees of freedom. In the 
massless case, only two of them interact with other parti- 
cles, namely the helicity-minus neutrino and the helicity- 
plus antineutrino. The other two degrees of freedom are 
sterile and have no impact at all. If neutrinos have a 
mass, their helicity eigenstates are no longer eigenstates 
of chirality and thus are no longer completely sterile. The 
strength of the weak interaction for "wrong-hclicity" neu- 
trinos is suppressed by a factor of about (m VT /2E VT ) 2 rel- 
ative to the usual rates, provided that m„ <C E v . In this 
case the quasi-sterile states remain essentially unexcited, 
leaving us with two degrees of freedom which are relevant 
for the cosmic expansion (for a recent discussion of cosmic 
abundances of right-handed neutrinos with masses up to 
200 keV see Enqvist et al. 1997). Conversely, for nonrela- 
tivistic neutrinos all four degrees of freedom interact with 
the full strength. Therefore, if the neutrino mass exceeds 
about 1 MeV so that they are nonrelativistic at the time 
of freeze-out, all four degrees of freedom would be fully 
excited and would thus contribute to the expansion rate 
of the universe 

A correct treatment of the intermediate-mass regime 
requires a detailed kinetic treatment of neutrino hclicities 
(e.g. Fields et al. 1997). Because this is a very significant 
complication we have limited ourselves to an approxima- 
tion where it is assumed that all available helicity states 
are always fully populated, an assumption that we shall 
refer to as the "helicity approximation" . It implies that the 



number and energy density of light (m„<l MeV) massive 
Dirac neutrinos is overestimated by a factor of 2. In the 
mass range from 1-10 MeV the number and energy den- 
sity is still overestimated, but now by a factor less than 
2; for mj,>10MeV our approximation yields the correct 
result. 

Because we have assumed kinetic equilibrium, the en- 
ergy density of the r neutrinos can be calculated from 
their number density. As an example we show in Fig. 3 
the evolution of the number and energy density, respec- 
tively, for a 10 MeV Majorana r neutrino. Even though the 
neutrino mass causes the universe to become quickly mat- 
ter dominated, we continue to parameterize the impact on 
the expansion rate by g* which thus no longer represents 
the effective number of relativistic degrees of freedom. It 
is merely a time-dependent proportionality factor between 
T 4 and the energy density according to the definition in 
Eq. (2). 

4 BBN with massive r neutrinos 

4-1 The number of effective light neutrino species 

It is common practice in the literature to parameterize 
the influence on BBN of a new particle by the equivalent 
number N e g of light neutrino species. This is a useful pa- 
rameterization for a particle which is relativistic during 
the entire period of BBN, i.e. down to T w 40 keV or un- 
til about 200 seconds after the big bang. Such a particle 
contributes a fixed amount to the parameter g* during the 
relevant epoch. On the other hand, a particle which is or 
becomes nonrelativistic during the BBN epoch contributes 
to the expansion rate in the way discussed in the previous 
section and displayed in Fig 2, i.e. its contribution to g* 
varies with time. Evidently, the effect of a massive r neu- 
trino cannot be mapped on a constant shift of or on a 
fixed modification of N c g. 

One can still define a useful N e g parameter indirectly 
by comparing the modification of the calculated abun- 
dance of a certain light element, for example 4 He, caused 
by a massive neutrino with that caused by a fixed mod- 
ification of g*. The synthesis of the light elements peaks 
at different times for the different elements. Therefore, 
defining N c g in this way depends on the chosen isotope. 
In Fig. 4 we show N e g based on the equivalent impact on 
the abundances of deuterium and the two helium isotopes. 
(Note that in the Dirac case the low-mass range does not 
converge to the correct value of 3 light families, but to 4. 
This represents the abovementioned helicity approxima- 
tion.) 

In the previous literature N c g has always been de- 
fined by the equivalent impact on 4 He. However, if one 
wishes to perform a more direct comparison with differ- 
ent sets of observational data, there is no reason to in- 
troduce N e g . We will present all of our results directly as 
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Fig. 4. Number of effective light neutrino families N e fi equivalent to the contribution of a massive neutrino. The value of N e g 
depends on the light element it is based upon. For Dirac neutrinos, our "helicity approximation" causes iV e ff to be overestimated 
for ra„ T < 10 MeV; this region has been hatched as a warning. 



a function of vh Vt in a form analogous to that used by 
Kawasaki et al. (1994). 

4-2 Validity of the approximated Boltzmann equation 

The N c s parameter remains very useful to compare our 
results with those of other authors who have dropped 
some or all of the simplifying assumptions made in Sect. 3 
to calculate the neutrino freeze-out. Fields et al. (1997) 
used the integrated version of the Boltzmann equation, 
but with correct Fermi-Dirac statistics for the neutrinos. 
Furthermore, they did not assume thermodynamic equi- 
librium for the electron neutrinos, but solved the coupled 
Boltzmann equation for all three neutrino families and 
took into account the influence of the increased v c num- 
ber density on the weak reaction rates. This changes the 
maximum value of N eS by about 0.5, but the limits on 
the allowed r neutrino mass are hardly affected, because 
below N c ff — 4 the results practically coincide. 

Hannestad & Madsen (1996a, b) solved the full Boltz- 
mann equation for all neutrino species without any sim- 
plifying assumptions. In addition they calculated the in- 
fluence of the deviation from equilibrium of the electron 
neutrino distribution function on the weak reaction rates. 
Their results deviate at most by 0.2 neutrino families com- 
pared to Fields et al. (1997). 

Dolgov et al. (1996) solved the integrated Boltzmann 
equation with the method of the pseudo chemical poten- 
tial only for the r neutrinos, but calculated the devia- 



tion from equilibrium of the electron neutrino distribu- 
tion function on the weak reaction rates. In this paper 
they do not present values for N c g 7 but only for AN c g 
relative to their earlier paper (Dolgov & Rothstein 1993) 
for the results without ^ c -heating. There, very large 
values for N c g were shown with A^™ c ax 11 and 

7V eff?max na ~ 8 - Their results for AN cS are opposite to 
those of the other groups. They argue that this is because 
Fields et al. (1997) did not take into account the devia- 
tion from equilibrium for the distribution function of the 
electron neutrinos, but Kainulainen (1996) showed in an 
update to Fields et al. (1997) that this would change their 
bounds by less than half a neutrino family. 

In summary, the modifications of N c g brought about 
by a more complete kinetic treatment of neutrino freeze- 
out and the weak interaction rates appear to be small 
enough to justify our simple treatment with the integrated 
Boltzmann equation. 

In the study of Kawasaki et al. (1994) the -/V e ff param- 
eter was not introduced, and they performed a rather com- 
plete numerical treatment of the v T kinetics. We have con- 
structed a contour plot with the observational data used 
by Kawasaki et al. (1994) in order to be able to compare 
directly with their Figure 5(e) which refers to the stable 
Majorana neutrino case. Our contours agree with theirs 
better than we would have expected, corroborating that 
our simple approximation is entirely justified for the task 
at hand. 
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Fig. 5. Isoabundance contours for the light elements if the neutrino mass is taken to be of Majorana type. In each panel, the 
contours have been drawn for a specific combination of observational data which are denoted according to the labels given in 
Table 1. (Left column: low 4 He. Right column: high 4 He. Top row: deuterium from ISM. Middle row: high deuterium from QAS. 
Bottom row: low deuterium from QAS.) The concordance regions for n and m Vr are shaded. 



4-3 Concordance Regions 

In Figs. 5 and 6 we present our results as isoabundance 
contours in the plane of baryon density vs. neutrino mass. 
In each panel, contours are drawn for 4 He (solid), deu- 
terium (dashed) and 7 Li (dotted). The contour levels cor- 
respond to the adopted observational ranges listed in Ta- 
ble 1. Each panel corresponds to a specific combination 
of observational results as described in the figure. Con- 
cordance regions in the r\-vn VT -plane are shaded. In the 



Dirac case (Fig. 6) the area below about lOMeV has been 
hatched because there our results are not quantitatively 
meaningful because of our helicity approximation. Qual- 
itatively, the Dirac and Majorana cases, of course, yield 
rather similar results. 

The Dirac r neutrino results have to converge to the 
ones for the Majorana case in the small mass regime 
(m UT < 1 MeV), since the "wrong-hclicity" states are not 
populated there. In the intermediate mass range from 1 - 
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Fig. 6. Isoabundance contours as in Fig. 5 if the neutrino mass is of Dirac nature. In the hatched area the results are not 
quantitatively correct because of our "helicity approximation" . 



10 MeV our approximation still overestimates the contri- 
bution of the Dirac r neutrino energy density to the cos- 
mic expansion rate; the correct results for the light clement 
abundances lie somewhere between the two cases . (For the 
He abundance in the presence of a massive Dirac r neu- 
trino with the correct treatment of the helicity states see 
Fields et al. 1997.) The correct treatment of both helicity 
states does not open up new concordance regions in this 
mass range. 

The most recent experimental limit on the r neu- 
trino mass is m„ r < 18 McV (ALEPH-collaboration, 



Passalacqua 1996), somewhat lower than the previously 
published limit 
(ALEPH-collaboration, Buskulic et al. 1995). Therefore, 
the concordance regions above this mass limit are already 
ruled out. 

If the classic ISM determination of the primordial deu- 
terium abundance were correct after all (top row in Figs. 5 
and 6) we would be left with a concordance region in the 
low-mass regime; in fact the neutrino mass may well be 
zero in this case. Even though the concordance region does 
not show up in the Dirac case because of our helicity ap- 
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proximation, it would be there because for low-mass neu- 
trinos the Dirac and Majorana cases are equivalent. The 
concordance region is almost independent of which 4 He 
result is used because of the constraints provided by the 
deuterium and lithium contours. A mass range between a 
few hundred keV and few tens of MeV is excluded. 

If the high QAS deuterium abundance is correct (mid- 
dle row in Figs. 5 and 6), the result depends on the as- 
sumed 4 He abundance. If the low value is right, the sit- 
uation is similar to that discussed in the previous para- 
graph, except that the inferred cosmic baryon abundance 
is smaller by a factor of 3. If the high helium abundance 
is the true value, then in the Majorana case there remains 
a small concordance region which actually requires a neu- 
trino mass somewhat below 1 MeV. In the Dirac case, the 
situation is unclear because of our helicity approximation. 

If the low QAS deuterium abundance is taken to be 
the primordial one (bottom row in Figs. 5 and 6), there 
is essentially no concordance region whatsoever, indepen- 
dently of the neutrino mass. This situation would consti- 
tute a true "crisis for BBN" as there would be no consis- 
tency between observations and calculations. Some other 
novel ingredient would be needed besides a massive neu- 
trino. 



5 Summary 

We have studied the influence of a massive r neutrino 
in the MeV range on the primordial abundances of the 
light elements. Several recent studies indicate that our 
simple approximate treatment of the Boltzmann collision 
equation for the neutrinos is quite adequate. We have 
focussed on a comparison between the calculated light- 
clement abundances with current observational data by 
virtue of contour plots in the t]-m VT -plane. This method 
avoids the use of an effective number of neutrino families 
N e g which cannot be defined simultaneously for all light 
element abundances. 

We find that the existence and location of concordance 
regions where all observations can be reproduced depends 
dramatically on the assumption which of the abundances 
currently offered as the "observed primordial" ones are 
actually correct. Depending on ones choice or judgment 
one may find that BBN is in crisis (no concordance region 
at all), that a nonvanishing neutrino mass is required, or 
that there are concordance regions even for a vanishing 
neutrino mass. 

In each case one finds that a neutrino mass between 
the experimental limit of 18 MeV and about 1 MeV (and 
a bit below) is ruled out. In this regard our conclusions 
agree with those of previous authors. In spite of the great 
uncertainty of the observational situation, this conclusion 
is rather stable because of the opposite curvature of the 
deuterium and the helium contours which peak in opposite 
^-directions for neutrino masses of a few MeV. 



Still, given the unclear observational situation one can- 
not be confident in the reliability of such limits. After all, 
if there were no concordance region whatsoever so that a 
new ingredient for BBN would be required, the possible 
impact of this new physics on the present results is entirely 
unknown. One clearly has to wait until the mutual incon- 
sistencies of the observationally inferred abundances have 
been sorted out before one can arrive at far-reaching cos- 
mological conclusions about the properties of elementary 
particles. 
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A Annihilation rates for massive r neutrinos 

The reaction rate for r neutrinos annihilating to fermions 
is given by the thermal average over the cross sec- 
tion a times the rclativistic invariant relative velocity v 
(Gondolo & Gelmini 1991) 



(<TV) 



4m 4 T TK\ (m„ T / T) 

f ds nvE VT E Pt Js-4ml T K 1 {y/s/T) . 

J Ami V 



(Af) 



Ki is the modified Bessel function of the second kind of 
order i as defined in Abramowitz & Stegun (1972), Ei are 
the particle energies, s = (P Vr +Pp T ) 2 is the center of mass 
energy, and P are the four momenta. We borrowed the ex- 
pression avE„ T Ep T from Dolgov & Rothstein (1993). The 
expressions given there are misprinted in several ways as 
confirmed by the authors (private communication). The 
correct expressions must read for Dirac neutrinos 
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and for Majorana neutrinos 
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Here, Pi are the three momenta, Gp the Fermi-constant, 
m a, Cv, a , and C\. a are the mass and the vector and ax- 
ial couplings for the final state fermions. Further we have 
defined w = (s 2 — ism 2 ,) 1 / 2 . The squared matrix element 
\M UT jy T ^ a a\ 2 is summed over the kinematically allowed re- 
action channels {v T v T — > e+e~, ^ e t> e , v^Dy) and averaged 
over the spins. 
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